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We discuss the calculation of resonant inelastic x-ray scattering (RIXS) with Xclaim, a multi- 
plaform code for the calculation of core-hole spectroscopy based on a ligand-field multiplet model. 
RIXS can be calcuated with a ligand-field model with ligands for the L 2 , 3 , M 2 , 3 , and Kp edges 
for d-valence metals. In addition, the program brings the possibility of fully diagonalizing the 
hamiltonian to classify the excited states. The command line use of xclaim is described for scripting 
purposes to help in automatizing the results that can be obtained with the graphical interface. 


Resonant inelastic x-ray scattering (RIXS) is an exper¬ 
imental technique for observing different types of excita¬ 
tions (dd, charge transfer, magnons, etc. ) that has at¬ 
tracted a lot of interest. ULi Amongst its advantages can 
be cited chemical and bulk selectivity, and the possibility 
to be used at high pressure. However, multiplet effects 
can make its interpretation not straightforward. A many- 
body model is necessary many times to take into account 
the character of the electronic structure of strongly cor¬ 
related materials, and the interaction with the core hole 
in the final states. 

Xclaim Hi is a code for the calculation of core-hole 
spectroscopy using a ligand-field multiplet model [El, [6| . 
The hamiltonian matrices and spectroscopy is calculated 
by a compiled fortran code with the input parameters be¬ 
ing set in a graphical interface in python (Fig. [T]). In this 
paper we descibe the extension of the code for calculating 
RIXS spectra and classifying the excited states. We also 
discuss the use of xclaim as a command line tool, that 
allows to reproduce previous calculations and automatize 
the process of calculating spectra. 

1. CALCULATION OF RIXS SPECTRA 

In the RIXS technique incoming x-ray photons {hw) 
are tuned in the vicinity of an absorption edge of the 
resonant ion and one measures the scattered intensity 
at for outgoing photons of energy fiuj'. The RIXS 
cross-section Q can be calculated through the Kramers- 
Heisenberg formula 
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1^) corresponds to the ground state of the initial state 
hamiltonian (without core-hole) and |/) labels the final 
states. The intermediate denominator corresponds to the 
Green’s function of the intermediate state hamiltonian 
Hi with F is the intermediate core-hole lifetime broaden¬ 
ing and the final 5 function imposes energy conservation. 
T and T' are the transition operators for the incoming 
and outgoing light, T = er and T' = e'r in the dipolar 



FIG. I: Main input window for xclaim. Atomic, crystal field, 
and charge transfer parameters for the calculation can be set 
in the input boxes. Pop-up windows can be used for the 
detailed modification of Coulomb and spin-orbit parameters, 
and the use of a general crystal-field in terms of Wybourne 
parameters. 

case for incoming and outgoing beams with polarization 
vectors e and e' 

A common scattering geometry is to use 90*^ between 
incident and scatterd beam to minimize the elastic scat¬ 
tering, and to use a and tt polarization for the incoming 
light. No polarization analysis is normally used for the 
outgoing beam, since it would atenuate too much the 
beam. 

The RIXS calculation done by full diagonalizations of 
the hamiltonian matrices can be subdivided into sev¬ 
eral independent processes. Each diagonalization can 
be separated into two independent processes for the up¬ 
per and lower eigenvalues (in the case of and 
RIXS three hamiltonian matrices need to be diagonal¬ 
ized). The lorentzian and gaussian convolutions of the 
RIXS map are also paralellized. The usage of the mul¬ 
tiprocessing python library allows to launch several in¬ 
dependent python processes from the python interpreter 
and doing the calculation in a computer with several pro¬ 
cessor cores will notably speed up the RIXS calculation. 
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FIG. 2: Output contour plot for a RIXS calculation for Co^^ 
at the Ka edge based on a crystal field model . 



2. EXCITED STATES 

In addition to getting information from the ground 
state it is possible to also characterize the excited states 
by doing an exact diagonalization of the hamiltonian. By 
selecting diagonalize in the spectroscopy tab the program 
will fully diagonalize the hamiltonian and show a text 
window with the energy positioning and information for 
the excited states (Fig. [3]). The text window shows for 
each quantum state the energy positioning the number 
of valence electrons, x, y and z components of the orbital 
angular momentum (L) and spin (S), the expectation 
value of the spin orbit operator for each state (labeled 
by L.S), the total spin (S^) = S{S + 1) (this is related 
to the spin multiplicity, i.e. singlet, doublet...), and the 
occupation of the different d orbitals: (ia, 2 _^ 2 , dxy, 

dyz^ and dzx- In the case of a calculation with ligands 
it will output the occupation of the valence and ligand 
shells. 

Information on excitations can be used in classifying 
the features of RIXS spectra as different kind of dd ex¬ 
citations or in the interpretation of pump-probe exper¬ 
iments, where the system undergoes a cascading decay 
through excited quantum states. Another inter¬ 

est for studying excitations comes from finding systems 
where crystal-field exctiations very close to the ground 
state allowing an easy switch the ground state by small 
changes in the ligand environment. 


3. COMMAND LINE INTERFACE 

In addition to using the graphical interface, it is pos¬ 
sible to run the program from the command-line and to 
repeat the calculation with the same parameters as set 
previously on the graphical interface, or modifying them. 
The file out/input.txt contains the parameters saved from 


FIG. 3: Window showing the results of diagonalizing the 
hamiltonian. 


the last execution of the program. It is worth noting that 
out/input.txt is rewritten each time the Calculate button 
is pressed or when the command line interface is run. The 
input file is saved as a python dictionary, between braces 
and different lines formated as ”labeF: value, pairs. Val¬ 
ues can be numerical or a string given between quotes. 
The out directory also contains other files depending on 
the type of calculation that is run: out/GS-info, txt (cal¬ 
culated expectation values for the ground state) and the 
files out/poles-*.txt (calulated spectra for different polar¬ 
izations). 

To use the comand line interface just type cp 
out/input.txt inputl.txt ; ./xclaim-gui.py -i inputl.txt 
The resulting spectra and ground state expectation val¬ 
ues will be in the directory out (its contents will be over¬ 
written) . A complete list of options can be seen by typing 
./xclaim-gui.py - -help The xclaim webpage Q contains 
more detailed information about the command line and 
rebroadening spectra. 
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